Superimposed Oscillations in the WMAP Data? 
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The possibility that the cosmic variance outliers present in the recently released WMAP multipole 
moments are due to oscillations in the primordial power spectrum is investigated. Since the most 
important contribution to the WMAP likelihood originates from the outliers at relatively small 
angular scale (around the first Doppler peak), special attention is paid to these in contrast with 
previous studies on the subject which have concentrated on the large scales outliers only (i.e. the 
quadrupole and octupole). As a physically motivated example, the case where the oscillations are 
of trans-Planckian origin is considered. It is shown that the presence of the oscillations causes 
an important drop in the WMAP % 2 °f about fifteen. The F-test reveals that such a drop has 
a probability less than 0.06% to occur by chance and can therefore be considered as statistically 
significant. 

PACS numbers: 98.80.Cq, 98.70.Vc 



I. INTRODUCTION 

The recently released WMAP data ,1] have confirmed 
the standard paradigm of adiabatic scale invariant pri- 
mordial fluctuations 0, 0. This paradigm can be jus- 
tified in the framework of inflation and can explain the 
most important cosmological observations 0, This 
remarkable success has led the cosmologists to take an 
interest in more subtle features of the WMAP multi- 
pole moments. In particular, recently, many studies have 
been devoted to the so-called cosmic variance outliers, i.e. 
points which lie outside the one sigma cosmic variance er- 
ror [S[l|l|l[l3[lJ[ll[lil. Th<3Se outliers have been 
considered as intriguing since the probability for their 
presence would be very small ^3- So far, all the studies 
have concentrated on the seeming lack of power at large 
scales, i.e. on the quadrupole and the octupole outliers. 
In the literature, two possibilities have been envisaged. 
In Ref. [HI , it has been argued that the outliers are not a 
problem at all, the crucial point being the way the proba- 
bility of their presence is estimated. In Refs. fA IsL M. ITof . 
it has been envisaged that the outliers could be a signa- 
ture of new physics even if it has also been recognized in 
these articles that the cosmic variance could be respon- 
sible for their presence. In particular, in Ref. Q, it has 
been proposed that the inflationary scale invariant initial 
power spectrum could be modified by some new physics 
such that a sharp cut-off at large scales appears while it 
remains unchanged elsewhere. It has been shown that 
this can cause a decrease of the \ 2 of order A% 2 ~ 2 — 4 
for one additional free parameter given by the scale of 
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the cut-off. 

However, as revealed by the Fig. 4 of Ref. [3, the 
main contribution to the WMAP % 2 does not come from 
the large scales but rather from scales which correspond 
to the first and second Doppler peaks (more precisely, 
according to Ref. |l4j , the three main contributions come 
from the angular scales £ ~ 120,200 and 340). In other 
words, if the presence of outliers is taken seriously into 
account, modifications of the standard power spectrum 
seem to be required on a different range of scales and, in 
any case, not only at very large scales. In addition, the 
small scale outliers can be above or below the theoretical 
error bar and, therefore, the required modifications do 
not seem to have the form of a systematic lack or excess 
of power. This naturally leads to the idea that the power 
spectrum could possess superimposed oscillations. The 
aim of this article is to study whether this idea has any 
statistical support. Of course, a physical justification 
for the presence of oscillations in the power spectrum is 
needed. Interestingly enough, such a justification exists 
precisely in the context of the theory of inflation flM IT5 
fl7| and we now discuss this question in more details. 

One of the main advantage of the inflationary scenario 
is that it permits to fix sensible initial conditions. Be- 
cause the Hubble radius was constant during inflation, 
the wavelength of a mode of astrophysical interest to- 
day was much smaller than the Hubble scale at the be- 
ginning of inflation whereas, without a phase of infla- 
tion, the same mode would have been a super-Hubble 
mode. Contrary to the super-Hubble case, the vacuum 
is defined without ambiguity in the sub-Hubble regime. 
This state is the starting point of the subsequent cosmo- 
logical perturbations evolution. This leads to a nearly 
scale-invariant spectrum for density perturbations, a pre- 
diction which is now confirmed to a high level of accu- 
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However, this remarkable success carries in itself a po- 
tential problem. In a typical model of inflation, the 
modes are initially not only sub-Hubble but also sub- 
Planckian, that is to say their wavelength is smaller than 
the Planck length m" 1 |l6lll7| . In this regime, the physi- 
cal principles which underlie the calculations of the power 
spectrum are likely not to be valid anymore. This prob- 
lem is specific to the perturbative approach and does not 
affect the background model. Indeed, since the energy 
density of the inflaton field is nearly constant during in- 
flation, we face in fact a situation where the wavelength 
of the Fourier modes is smaller than the Planck length 
whereas the energy density which drives inflation can still 
be well below the Planck energy density. The problem de- 
scribed above mainly concerns the modes of cosmological 
(still in the linear regime) interest today which are, in the 
inflationary paradigm, a pure relic of the trans-Planckian 
regime. It should be added that the scale at which the 
new physics shows up is not necessarily the Planck length 
but could be everywhere between this scale and the Hub- 
ble radius. In this article we denote the new energy scale 
M c and assume that a = H/M c is a free parameter 
{H denotes the Hubble parameter during inflation). If 
the physics is different beyond the scale M c , this should 
leave some imprints on the spectrum of inflationary cos- 
mological perturbations and therefore modify the Cosmic 
Microwave Background (CMB) anisotropics |15lll6lll7| . 

The next question is how to calculate these modifica- 
tions. Roughly speaking, the calculation of the power 
spectrum of the fluctuations reduces to the calculation 
of the evolution of a free quantum scalar field in a time- 
dependent background. Various methods have been used 
to model the physics beyond M c like changing the free 
dispersion relation Il6t Il7l Il8l Il9l. using stringy uncer- 
tainty relations poll2lLl22ll23ll24 l25l| or noncommutative 
geometry [2(| . Yet another approach has been to assume 
that the Fourier modes are "created" when their wave- 
length equals the critical scale M" 1 [2 ^ j2p j 2 p |30ll3l | . A 
generic prediction, first made in Refs. |16u17|. which ap- 
pears to be independent of the settings used to model the 
new physics, is the presence of oscillations in the power 
spectrum (of course, the detailed properties of these os- 
cillations do depend on the model utilized). This can 
be understood easily. In the standard calculation, ini- 
tially, the scalar field is just given by an in-going wave. 
If the evolution proceeds according to the WKB approxi- 
mation [32, at Hubble scale exit the scalar field will just 
differ by a phase which will drop out when the power 
spectrum is calculated. On the other hand, if the WKB 
approximation is violated at some time before Hubble 
scale exit or if the modes are created at a fixed length in 
some a-vacuum, then the scalar field will be a combina- 
tion of on-going and out-going waves. This combination 
gives the oscillations in the power spectrum. To put it 
differently, one can say that the existence of a preferred 
scale plus the standard inflationary theory generically im- 
ply a power spectrum given by a nearly scale invariant 
component plus superimposed oscillations. 



The aim of this article is to use the previous generic 
trans-Planckian prediction as a tool to analyze whether 
the WMAP multipole moments exhibit oscillations orig- 
inating from oscillations in the primordial spectrum. It 
should be stressed again that, in this work, the trans- 
Planckian effects are only considered as an example of 
models where oscillations could show up. Indeed, the 
analysis presented in this article would still be useful for 
any alternative model which predicts a similar oscillatory 
pattern. 

This paper is organized as follows. In the next section, 
we briefly recall the standard calculation of the inflation- 
ary spectrum when the trans-Planckian effects are taken 
into account. In Sect. IIIII we compare the theoretical 
predictions with the WMAP data. Finally, we discuss 
the data analysis and give our conclusions in Sect. HVl 



II. THE TRANS-PLANCKIAN POWER 
SPECTRA 

In this article, we assume that inflation of the spa- 
tially flat FLRW spacetime is driven by a single scalar 
field 4'o( r l)- Scalar perturbations of the geometry can 
be characterized by the gauge-invariant Bardeen poten- 
tials $ and fluctuations in the scalar field are character- 
ized by the gauge-invariant quantity 8<p. If <f>' 7^ (a 
prime denotes derivative with respect to conformal time 
77), then everything can be reduced to the study of a 
single gauge-invariant variable (the so-called Mukhanov- 
Sasaki variable) defined by 33j v = a(54> + cj)' ( $/Tl), 
where TL = aH is the conformal Hubble parameter. 
In fact, it turns out to be more convenient to work 
with the variable /i s defined by ^ s = —\^2kv, where 
k = Sir/m 2 ^ Density perturbations are often charac- 
terized by the so-called conserved quantity £ [34| defined 
by C = (2/3) (TlT 1 *' + $)/(l + w) + $, where w is the 
equation of state parameter. The quantity fi s is related 
to C by fi s = -2a^yQ, where 7=1- H'/H 2 . On the 
other hand, the primordial gravitational waves are de- 
scribed by the transverse and traceless tensor h^, the 
Fourier transform of which is /i T (up to a polarization 
tensor) . Both types of perturbations obey the same type 
of equation of motion, namely the equation of a paramet- 
ric oscillator 



(1) 



with lu 2 s = k 2 - (a^7)"/(a % /7), u\ = k 2 - a" /a, and 
k is the wavenumber of a given Fourier mode. Fi- 
nally, the quantities of interest for computing the CMB 
anisotropies are the power spectra which read [34[ 
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In order to compute k P^(k) and fc 3 P/j(fc), one must in- 
tegrate the equation of motion Q and specify what the 
initial conditions are. The integration of the equations of 
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motion is possible for a large class of inflationary models 
provided they satisfy the slow-roll conditions j35'] . The 
slow-roll approximation is controlled by a set of parame- 
ters given by e = 7, S = /(H(/>' )-l and ^= (e'-6')/H. 
The slow-roll conditions are satisfied if e and S are much 
smaller than one and if £ = 0(e 2 , S 2 , eS). At first order, 
the parameters e and 6 can be considered as constant. 

Here, the initial conditions are fixed under the assump- 
tion that the Fourier modes never penetrate the trans- 
Planckian region. In other words, a Fourier mode is sup- 
posed to "appear" when its wavelength becomes equal 
to a new fundamental characteristic scale £ c . The time 
r/fe of mode "appearance" with comoving wavenumber fc, 
can be computed from the condition 



= 4 = ^, 



(3) 



which implies that rjk is a function of k. This has 
to be compared with the standard inflationary calcula- 
tions where the initial time is taken to be rjk = —00 
for any Fourier mode k and where, in a certain sense, 
the initial time does not depend on k (see, however, 
Ref. |36(). In the framework of trans-Planckian infla- 
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tion, a crucial question is in which state the Fourier 
mode is created at the time rjk- Here, we consider 
the most general conditions,also called a truncated a- 
-vacuum EHM IS E3, El El 
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The coefficients and are a priori two arbitrary com- 
plex numbers satisfying the condition \ck\ 2 — \dk\ 2 = 1. 
Since there are two energy scales in the problem, namely 
the Hubble parameter H during inflation and the new 
scale M c , the final result will be expressed in terms of 
the ratio a g = H/M c , which is a small parameter. As a 
result, we typically expect that Cfc = 1 + ya + ■ ■ ■ and 
dk = xa + • • • since the initial conditions are expressed 
at time r\k |3Cj. The parameters x and y are considered 
as free parameters that are not fixed by any existing well- 
established theories (see also Ref. except, of course, 
that they should be such that the relation |cfe| 2 — |dfc| 2 = I 
is satisfied. One easily shows that this implies y + y* = 
at leading order in . Expanding everything in terms of 
<7„ , one arrives at l31l 
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where ip is the argument of the complex number x, i.e 
x = \x\e tv . The constant C is given by C = 7 E + In 2 — 2, 
7 E being the Euler constant. The scales fc* is the pivot 
scale Ej, . The parameter <r g and the scale factor a 
are evaluated at the time rj during inflation which is a 
priori arbitrary but, and this is the important point, does 
not depend on k. In the following we will choose this time 
such that fc*/a = M c . One sees that we obtain a scale- 
invariant spectrum plus logarithmic corrections in k the 
amplitude of which is determined by the slow-roll param- 
eters. In addition, superimposed oscillations coming from 
the trans-Planckian initial conditions of Eqs. J3J and © 
appear. The magnitude of the trans-Planckian correc- 
tions are linear in the parameter a and their amplitude 
is given by | a; j cr , a result in agreement with Ref. |30| . 
Let us also remark that the parameter y does not appear 
at the leading order considered in Eqs. I© and J7J). The 



wavelength of the oscillations can be expressed as 
k e 

in agreement with Ref. [43. The presence of the factor 
\x\ in the amplitude of the trans-Planckian corrections 
plays an important role in what follows. As stressed in 
Ref. [H, Refs. HHHEI assume that \x\ = 1. This has 
the consequence that high frequency waves have neces- 
sarily a small amplitude [see Eq. JSJl]. The factor |a;| al- 
lows to break this degeneracy. As a result the amplitude 
and the wavelength of the corrections become two inde- 
pendent quantities. The data analysis can then proceed 
in a larger parameter space. This explains why we can 
obtain results different from those derived in Ref. E|J 
without being in contradiction with that article. How- 
ever, let us also remark that a coefficient \x\ 7^ 1 implies 
that the corresponding initial state cannot belong to the 
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class of vacua often considered in previous works as, for 
instance, the instantaneous Minkowski vacuum state or 
the initial state considered in Ref. [23]. In this case, and 
as it will be discussed below, we face the issue of back- 
reaction [371 l38l l39i l40j . Finally, let us also notice that 
the power spectra depend on both parameters e and S 
and not only on the parameter e, contrary to some re- 
sults recently obtained in the literature. 

The initial power spectra of Eqs. iJBJ and Q are related 
to the CMB anisotropy through the multipole moments 
which are in turn defined through the two-point correla- 
tion function of the temperature fluctuations. Explicitly, 
we have 

/^( ei )f (e 2 )\ = £ *±lc t P t (coB6) , (9) 

where 9 is the angle between the two directions ei and e2. 
The WMAP satellite has also measured the polarization 
and therefore we will also be interested in the "TE" cross 
polarization multipole moments CJ E 01 ■ The (temper- 
ature) multipole moments are the sum of the scalar con- 
tribution and of the tensor contribution, Ci = Cf + CJ. 
From Eqs. J^jJ and Q, we see that in the context of 
slow-roll inflation, with or without trans-Planckian cor- 
rections, the contribution of gravitational waves cannot 
be put to zero arbitrarily. Indeed, the ratio of the two 
types of perturbations is predicted by the form of the 
power spectra and reads k 3 Ph/k 3 P^ — 16e(l+- • • ), where 
the dots stand for the trans-Planckian corrections. The 
fact that the trans-Planckian physics modifies the con- 
sistency check of inflation was first noticed in Ref. [48| . 
Finally, we notice that, since we work at first order in the 
slow-roll parameters, there is no running of the spectral 
indices. 

The derivation of the trans-Planckian corrections in 
the power spectra © and 0) assumes that the back- 
reaction effects are not too important. For consistency, 
the energy density of the perturbations must be smaller 
or equal than that of the inflationary background. This 
leads to the condition |a;| < ^/3Trm pl /Al c , an estimate 
which is in agreement with that derived in Ref. |3dj| . In 
order to put numerics on the above constraint, we can use 
the large scale approximation of the multipole moments 

C|c|^ +00 ^(fcn ss )fc 3 Pcy> (10) 

where je is a spherical Bessel function of order £ and r\ ss 
is the distance to the last scattering surface. This ap- 
proximation is valid for £ < 10 and also requires a van- 
ishing cosmological constant with no integrated Sachs- 
Wolfe effect. Neglecting the logarithmic corrections in 
the power spectrum gives l{£ + 1)C| ~ 27J 2 /(25em 2 i ). 
On the other hand, since the COBE-WMAP normaliza- 
tion is 6T 2 C 2 S x 10 12 /(2tt) ~ 10 3 , where T ~ 2.72K, this 
implies that H ~ m pl y/el0~ 4 . Finally, using the fact 



that M c = H/dQ, one arrives at 

2 

\x\a < 10 4 x ^, (11) 

where, in order to derive an order of magnitude estimate, 
the unimportant factors of order one have been neglected. 
It is important to emphasize that the above constraint 
is only a sufficient condition, but by no means, unless 
proven otherwise, a necessary condition for the validity 
of the power spectra calculations. 

III. COMPARISON WITH THE CMB DATA 

In this section, the compatibility of the trans- 
Planckian models described in the previous section with 
the WMAP data is studied. More precisely, we now use 
the trans-Planckian power spectra derived previously to 
test the idea that the cosmic variance outliers could be 
due to oscillations in the WMAP data. In order to com- 
pute the multipole moments Cg and C 'J E . we have used 
a modified version of the CAMB code based on its 
slow- roll inflation module [Hof. Following the analysis 
of Ref. [8|, the parameters space that we consider is h, 
^bi ^dm, z ic (or t), e, 5, P S caiar and the trans-Planckian 
parameters a Q , x. The parameter h is the dimension- 
less Hubble constant, fib the amount of baryons, fidm 
the amount of dark matter and z re the redshift of reion- 
ization. P S caiar is the normalization of the scalar power 
spectrum, and as mentioned above, it also fixes the nor- 
malization of the gravitational waves. We have restricted 
ourselves to fiat models, i.e. the cosmological constant is 
given by = 1 — fidm — fib ■ The choice ip = has been 
made since we have checked that the value of if has no 
real influence for our purpose. Therefore, we deal with a 
9-dimensional parameter space. In order to estimate the 
parameters, we have used Monte-Carlo methods imple- 
mented in the CDSMDMC code j^U with our modified CAMB 
version, together with the likelihood code developed by 
the WMAP team 0,0. 

We have first checked that, without trans-Planckian 
effects, i.e. by setting \x\a a — in the expressions of the 
scalar and tensor power spectra, the standard results are 
recovered. The best WMAP fit given in Ref. (3 is well 
reproduced, as well as the values for the slow-roll param- 
eters derived in Refs. 0,0,E2- Concerning the likelihood 
computation, by setting the input parameters to the best 
fit values given in Ref. we obtain x 2 ~ 1430.92 for 
1342 degrees of freedom, also in good agreement with 
Ref. (see the first line in Tabled). 

In a second step, we have taken into account the oscilla- 
tions present in the primordial power spectrum. Roughly 
speaking, the multipole moments Cg arc given by the 
convolution of a spherical Bessel function with the initial 
power spectrum [see Eq. (I10f) ]. It is therefore not ob- 
vious that the oscillations in the primordial spectra will 
be transferred into oscillations in the multipole moments. 
As shown in Ref. [53| , fine structures in the initial power 
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0.001 
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-0.013 
0.061 


0.028 
-0.083 


0.97 
1.04 


1430.92/1342 
1415.38/1340 



TABLE I: Best fit parameters from the WMAP data for the standard model of inflation compared to the best fit parameters 
obtained for the trans-Planckian model. The pivot scale and the time r/ (see the text) are chosen such that k*Ja a = M c = 
O.OlMpc -1 . The quantities ei and £2 are two different slow-roll parameters, defined in Ref. and used in Ref. |f| to express 
the observational constraints obtained from the WMAP data. These parameters are linearly related with e and S. 



spectra are strongly broadened and damped by the trans- 
fer function. However, the oscillations considered here 
are not localized at a given momentum scale but are 
rather spread all over the A; -space. One may therefore 
expect constructive interferences to appear at some an- 
gular scale £ even in presence of the damping mentioned 



above. From the analytical point of view, the integral 
in Eq. I|10|l can be evaluated exactly for the case 2e = 5 
(i.e. n s — 1) and this will allows us to gain some intu- 
ition with regards to the behavior of the trans-Planckian 
multipolc moments. They can be expressed as 
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where T stands for the Euler's integral of second kind |55j . 
and with £ < 10. Although this equation is not specially 
illuminating, we can already conclude that, due to the 
presence of the Eulcr functions with complex arguments, 
the oscillations are indeed passed to the multipole mo- 
ments. This is already an important conclusion since 
this means that the oscillations in the primordial power 
spectrum are not killed by the convolution. In the limit 
e/<r 3> £, the above formula simplifies a lot and reads 



£{t + l)C\ 



2H 2 
25em 2 



■(l-2eK 1 + V¥ 
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1 + e In ■ 



a M c r is 
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(13) 

One sees that, at small £, these oscillations are damped 
by a factor of (e/<7 ) 5 / 2 , for a given primordial amplitude 
|ar|<r . Note that at larger I, this effect may be com- 
pensated by the factor £(£ + 1). Therefore, in practice, 
we expect the oscillations to appear only at relatively 
small scales and to be absent at large angular scales. Of 
course, at large £, or for not too small values of a , the 
above equation quickly becomes invalid and an accurate 
estimation can be made only with the help of numerical 
calculations. 



These computations are not trivial and require the 
CAMB accuracy to be increased in order to correctly trans- 
fer the superimposed oscillations in the power spectra to 
the multipole moments. Indeed, the CAMB default compu- 
tation target accuracy of 1% for the scalars is only valid 
for non-extreme models. Since we precisely consider fine 
structures in the power spectra, the corresponding mod- 
els do not enter into this class and modifications in the 
default values of the CAMB code parameters are therefore 
required for small values of a . Without these modifica- 
tions (i.e. with the default accuracy), random structures 
in the multipole moments were found to appear at large 
scales for models with a < 1CP 2 . The modified accuracy 
of the code has been chosen such that it allows to treat 
models with sufficiently small values of a (in practice, 
up to <j Q = 10~ 4 ) and to avoid prohibitive computation 
time. Then, the strong damping of the high frequency 
oscillations is recovered at large scale. In Fig. ^ the 
difference between the trans-Planckian and inflationary 
angular power spectra, CJ PL — C™ MAP , has been plot- 
ted for three different values of a . As mentioned above, 
the previous qualitative considerations are recovered: the 
smaller a and I, the more damped the oscillations in the 
angular power spectrum. For values of a < 1CP 3 , the 
remaining oscillations are even washed out at small £ and 
only show up at the rise of the first acoustic peak. We 
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FIG. 1: Difference between the trans-Planckian and the stan- 
dard inflationary multipole moments, plotted for three values 
of <t and for an amplitude of \x\a = 0.3. A strong damping 
of the oscillations is observed at large angular scales. 



also notice a similar behavior for CJ E (see Fig. EJ. 

Having checked that the numerical calculations are well 
under control, one can now move to parameters estima- 
tion. The main result of this article is that, with the 
oscillations taken into account, it is possible to decrease 
the x 2 significantly. The best fit that has been found 
by C0SM0MC and the result of its Markov chains explo- 
ration of the parameter space is summarized in Table [I] 
It leads to \ 2 — 1415.38 for 1340 degrees of freedom, 
i.e. Ax 2 — 15 compared to WMAP one. The corre- 
sponding multipole moments Cg and CJ E together with 
the WMAP data are represented in the bottom panels of 
Fig. El The reason for such an important improvement of 
the x 2 is clear from these figures: the presence of the os- 
cillations allows a better fit of the cosmic variance outliers 
at small scales. In order to make this statement more 
quantitative, the difference between the cumulative x 2 
with and without the oscillations (in this last case this is 
nothing but the WMAP cumulative x 2 ) has been plotted 
in Fig. |21 as function of the angular scale £. This permits 
a direct comparison with the Fig. 4 of Ref . 0] . Clearly, 
the significant decrease of the x 2 comes from the outliers 
around the first acoustic peak which are well fitted by the 
oscillating component of our initial power spectra. No- 
tice also that the outliers at large scales, £ < 50, can not 
be well fitted in this model due to the strong damping de- 
scribed above [see Eq. (fH51) ]. But, as already mentioned 
several times, this is not a problem since the outliers that 
contributes the most to the x 2 are n °t located at large 
scales but around the first acoustic peak. 

To conclude this section, let us remark that we have 
not given confidence values for the best fit model pa- 
rameters in Table H] This is due to the fact that the 



high-accuracy CAMB version, which is necessary to obtain 
physical results when the oscillations are included, takes 
a few minutes to compute only one model instead of a 
few seconds in the standard inflationary case. Although 
minima in the x 2 hypersurface can be "rapidly" found 
with C0SM0MC , the Markov chains take a much longer 
time to scan the "one-sigma" neighborhood. In addition, 
as studied in Ref. ^(|, the likelihood function looks like 
a hedgehog in the parameters space forbidding a smooth 
convergence of the Markov chains around a given x 2 min- 
imum. 



IV. DISCUSSION AND CONCLUSION 

Given the previous result, Ax 2 — 15, the first question 
that one may ask is whether the improvement is statisti- 
cally significant? We have introduced two new parame- 
ters, <7 and x, and one may wonder whether it was worth 
it, given the decrease observed in the x 2 - Indeed, it could 
just be due to statistical fluctuations in the larger param- 
eters space. Since it is not possible to directly compare 
the x 2 °f models with different degrees of freedom, i.e. 
with different number of parameters, we need another re- 
liable statistical test. The F-test |5Jj is an efficient tool 
to deal with this problem. This test works as follows. 
Given a first fit with n 1 degrees of freedom and with a 
x\ and a second fit with n 2 < n 1 degrees of freedom and 
a new x 2 < x\ > the F-test gives the probability that the 
decrease in the x 2 is only due to statistical fluctuations 
and not to the fact that the underlying model is actually 
a better fit. Generally, one considers that the improve- 
ment is significant if this probability is less than a few %. 
For our case, the F-probability is F pro b a — 0.06%. 

A possible loophole would be that, since the F-test is 
only valid for Gaussian statistics, it could not be applied 
to the case at hand since the multipole moments obey a 
different statistics. However, since our result comes from 
relatively small scales, the central limit theorem applies 
and the statistics should not deviate too much from the 
Gaussian one. Therefore, we expect the F-test to give 
a rather fair estimate. Another test has been to check 
that a spectrum with a different oscillatory pattern, e.g. 
with a fc-dependence instead of In k in the cosine function, 
does not provides us with a fit as good as the trans- 
Planckian one, in the F-test sense and for a same amount 
of computed models with C0SM0MC . This is indeed the 
case since we have found A% 2 ~ 4 only (with two new 
parameters) which is not very statistically significant. 

Let us now analyze and discuss the result itself in more 
details. The first remark is that the baryons contribution 
is standard since we have Slbh 2 ~ 0.025. It is therefore 
compatible with Big Bang Nucleosynthesis even if this 
is slightly less compatible than the value obtained with- 
out oscillations but this difference does not seem to be 
significant. In the same way, the value of fidm is not sig- 
nificantly modified compared to the standard one. The 
same is true for the energy scale of inflation, i.e. for 
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FIG. 2: Top left and right panels: the best fit obtained for trans-Planckian angular TT and TE power spectra compared to 
the best fit standard inflationary power spectra. Bottom left and right panels: the same best fit trans-Planckian TT and TE 
power spectra but this time compared to the binned WMAP data 



■Pscaiar- On the other hand, we see that the value of Q\ 
is particularly high, as the values of h and the optical 
depth are. However, this still seems to be admissible, 
in particular with the SNIa measurements. However, the 
slow-roll parameters differ significantly compared to their 
standard values |J>|5|,|52j and, as a consequence, the spec- 
tral index is modified. This is not surprising and should 
even be expected: since the trans-Planckian parameters 
seem to be as statistically significant as the slow- roll ones, 
their inclusion in the best fit search should have an effect 
on the determination of the slow-roll parameters. The 
fact that the best fit cosmological parameters with or 
without oscillations are not exactly the same illustrates 
the fact that the determination of these cosmological pa- 
rameters depends on the shape of the primordial power 
spectra |57j . 

Let us now study the trans-Planckian parameters in 



more details. First of all, we see that the best fit is 
such that the preferred scale M c is almost four orders 
of magnitude higher than the Hubble scale during in- 
flation. However, this does not provide us with a mea- 
sure of M c since the energy scale of inflation cannot be 
deduced directly from the data. Nevertheless, a con- 
straint on H has been given in Ref. [j| and leads to 
M c /m pl < 0.1. On the other hand, from Table U one 
obtains that 4 x lOV^/y^e ~ 10~ 2 . This means that 
the trans-Planckian amplitude, |a;|cr = 0.29, is not com- 
patible with the requirement of negligible backreaction 
effects [see Eq. I|llf) ]. This is certainly a major diffi- 
culty not for the presence of oscillations in the multipole 
moments (since we have proven that this is statistically 
significant anyway) but rather for the physical interpre- 
tation of those oscillations in terms of trans-Planckian 
physics. At this point, two remarks are in order. First, 
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FIG. 3: Cumulative residual Ax 2 {£) between the trans- 
Planckian fit and the WMAP best fit. The contributions 
stemming from the temperature and the polarization angular 
power spectra are also plotted. The final A\ 2 improvement 
is of order 15. 

the previous calculation does not predict what the mod- 
ifications coming from the inclusion of the backreaction 
effects are and, in particular, it does not tell that the 
backreaction effects will modify the power spectra. It 
just signals when the backreaction effects must be taken 
into account. Therefore, despite the fact that we are not 
able to prove it, it could very well be that the backreac- 
tion effects do not modify the fc-dependence of the oscil- 
lations (for instance, it could only renormalize the val- 
ues of the trans-Planckian parameters). Clearly, this is a 
very difficult technical question since a check of the above 
speculation would require a calculation at second order 
in the framework of the relativistic theory of cosmological 
perturbations with the trans-Planckian effects taken into 
account. This is beyond the scope of the present article. 
Secondly, the physical origin of the oscillations could be 
different. For instance, in Refs. |58l |59j, initial power 
spectra with oscillations were also found but with differ- 
ent physical justifications. In the model of Ref. |5q|. the 
presence of the oscillations is due to non-standard initial 



conditions in the framework of hybrid inflation, while 
in Ref. |59| a sudden transition during inflation is in- 
volved. Of course, it remains to be proven that a different 
model could produce an oscillatory pattern similar to the 
trans-Planckian one but since, on a purely phenomcno- 
logical level, it seems that the significant increase of the 
likelihood requires superimposed oscillation in the initial 
power spectra, this is a question that is certainly worth 
studying. Let us also signal that a primordial bouncing 
phase, as described for instance in Ref. |(3(j , may also do 
the job. Finally, let us remind that the main goal of the 
paper was to study whether the idea that oscillations are 
present in the WMAP data can be statistically demon- 
strated and that the trans-Planckian power spectra used 
before were just one possible example. 

As a conclusion, a word of caution is in order. From 
the previous considerations, it is clear that, in the ab- 
sence of outliers around the first acoustic peak, there 
would be no reason, in the F-test sense, to include oscil- 
lations with a large amplitude (i.e. |x| ^ 1) in the power 
spectra. Therefore, the future WMAP data release will 
be extremely important to check whether the presence 
of these outliers is confirmed or if they are just observa- 
tional artifacts. However, if the last possibility turns out 
to be true, it is clear that it would still be worth look- 
ing for small trans-Planckian effects (i.e. \x\ — 1) in the 
future very high-accuracy Planck data and the ideas put 
forward in the present article would still be useful for this 
purpose. 



Acknowledgments 

We wish to thank G. Hebrard, A. Lecavelier des 
Etangs, M. Lemoine, M. Peloso, P. Peter, S. Prunet and 
B. Revenu for helpful comments and/or careful reading of 
the manuscript. We are especially indebted to S. Leach 
for his help and many enlightening discussions. It is also 
a pleasure to thank R. Trotta for his useful advice on the 
CAMB and CDSMOMC armory. We would like to thank the 
CINES for providing us one year of CPU-time on their 
SGI-O3800 supercomputers. 



[1] C. L. Bennet et al. , Astrophys. J. Suppl. 148, 1 (2003), 

astro-ph/0302207. 
[2] G. Hinshaw et al. , Astrophys. J. Suppl. 148, 135 (2003), 

astro-ph/0302217. 
[3] L. Verde et al. , Astrophys. J. Suppl. 148, 195 (2003), 

astro-ph/0302218. 
[4] H. V. Peiris et al Astrophys. J. Suppl. 148, 213 (2003), 

astro-ph/0302225. 
[5] S. Leach and A. Liddle, Phys. Rev. D68, 123508 (2003), 

astro-ph/0306305. 



[6] G. Efstathiou, Mon. Not. Roy. Astron. Soc. 343, L95 

(2003), astro-ph/0303127. 
[7] C. Contaldi, M. Peloso, L. Kofman and A. Linde, JCAP 

0307, 005 (2003), astro-ph/0303636. 
[8] J. M. Cline, P. Crotty and J. Lesgourgues, JCAP 0309, 

010 (2003), astro-ph/0304558. 
[9] B. Feng and X. Zhang, Phys. Lett. B570, 145 (2003), 

astro-ph/0305020; B. Feng, M. Li, R.-J. Zhand and 

X. Zhang, Phys. Rev. D68, 103511 (2003), astro- 

ph/0302479. 



9 



[10] S. Tsujikawa, R. Maartens and R. H. Brandenberger, [36 
Phys. Lett. B574, 141 (2003), astro-ph/0308169. 
J.-P. Luminet et al. , Nature 425, 593 (2003). [37 
Q.-G. Huang and M. Li, JCAP 0311, 001 (2003), astro- 
ph/0308458; JHEP 06, 014 (2003), hep-th/0304203. [38 
G. Efstathiou, Mon. Not. Roy. Astron. Soc. 346, 
L26 (2003), astro-ph/0306431; 348, 885 (2004), astro- [39 
ph/0310207. 

D. N. Spergel et al. , Astrophys. J. Suppl. 148, 175 [40 
(2003), astro-ph/0302209. 

X. Wang, B. Feng and M. Li, astro-ph/0209242. [41 
J. Martin and R. H. Brandenberger, Phys. Rev. D63, 
123501 (2001), hep-th/0005209. [42 
R. H. Brandenberger and J. Martin, Mod. Phys. Lett. [43 
A16, 999 (2001), astro-ph/0005432. 

J. C. Niemeyer, Phys. Rev. D63, 123502 (2001), astro- [44 
ph/0005533. 

M. Lemoine, M. Lubo, J. Martin and J. P. Uzan, [45 
Phys. Rev. D65, 023510 (2002), hep-th/0109128. 
A. Kempf, Phys. Rev. D63, 083514 (2001), astro- [46 
ph/0009209. 

R. Easther, B. R. Greene, W. H. Kinney and G. Shiu, [47 
Phys. Rev. D64, 103502 (2001), hep-th/0104102. 
R. Easther, B. R. Greene, W. H. Kinney and G. Shiu, [48 
Phys. Rev. D67, 063508 (2003), hep-th/0110226. 
A. Kempf and J. C. Niemeyer, Phys. Rev. D64, 103501 [49 
(2001), astro-ph/0103225. 

R. H. Brandenberger and P. M. Ho, Phys. Rev. D66, [50 
023517 (2002), hep-th/0203119. [51 
S. F. Hassan and M. S. Sloth, Nucl. Phys. B674, 434 
(2003), hep-th/0204110. [52 

F. Lizzi, G. Mangano, G. Miele and M. Peloso, JHEP 
0206, 049 (2002), hep-th/0203099. [53 
U. H. Danielsson, Phys. Rev. D66, 023511 (2002), hep- 
th/0203198. [54 
R. Easther, B. R. Greene, W. H. Kinney and G. Shiu, 
Phys. Rev. D66, 023518 (2002), hep-th/0204129. [55 

G. L. Alberghi, R. Casadio and A. Tronconi, Phys. Lett. 
B579, 1 (2004), gr-qc/0303035. [56 
C. Armendariz-Picon and E. A. Lim, JCAP 12, 006 
(2003), hep-th/0303103. 

J. Martin and R. H. Brandenberger, Phys. Rev. D68, 
063513 (2003), hep-th/0305161. [57 
J. Martin and D. J. Schwarz, Phys. Rev. D67, 083512 
(2003), astro-ph/0210090. [58 
V. F. Mukhanov, H. A. Feldman and R. H. Branden- 
berger, Phys. Rept. 215, 203 (1992). [59 
J. Martin and D. J. Schwarz, Phys. Rev. D57, 3302 
(1997), astro-ph/9704049. [60 

E. D. Stewart and D. H. Lyth, Phys. Lett. B302, 171 
(1993). 



D. J. H. Chung, A. Notari and A. Riotto, JCAP 0310, 
012 (2003), hep-ph/0305074. 

T. Banks and L. Mannelli, Phys. Rev. D67, 065009 
(2003), hep-th/0209113. 

M. B. Einhorn and F. Larsen, Phys. Rev. D67, 024001 
(2003), hep-th/0209159. 

N. Kaloper et al. , JHEP 0211, 037 (2002), hep- 
th/0209231. 

K. Goldstein and D. A. Lowe, Nucl. Phys. B669, 325 
(2003), hep-th/0302050. 

H. Collins, R. Holman and M. R. Martin, Phys. Rev. 

D68, 124012 (2003), hep-th/0306028. 

H. Collins and R. Holman, hep-th/0312143. 

J. Martin and D. J. Schwarz, Phys. Rev. D62, 103520 

(2000), gr-qc/9704049. 

J. Martin, A. Riazuelo and D. J. Schwarz, Astrophys. J. 
543, L99 (2000), astro-ph/0006392. 

L. Bergstrom and U. H. Danielsson, JHEP 0212 038 

(2002) , hep-th/0211006. 

0. Elgar0y and S. Hannestad, Phys. Rev. D68, 123513 

(2003) , astro-ph/0307011. 

A. Kogut et al. , Astrophys. J. Suppl. 148, 161 (2003), 
astro-ph/0302213. 

L. Hui and W. Kinney, Phys. Rev. D65, 103507 (2002), 
astro-ph/0109107. 

A. Lewis, A. Challinor and A. Lasenby, Astrophys. J. 
538, 473 (2000), astro-ph/9911177, http://camb.info 



S. Leach, http : //astronomy . Sussex. ac . uk/~sleach/inf lation/ cam 
A. Lewis and S. Bridle, Phys. Rev. D66, 103511 (2002), 
astro-ph/0205436, http: //cosmologist . inf o/cosmomc 
V. Barger, H.-S. Lee and D. Marfatia, Phys. Lett. B565, 
33 (2003), hep-ph/0302150. 

W. Hu and T. Okamoto, Phys. Rev. D69, 043004 (2004), 
astro-ph/0308049. 

D. J. Schwarz, C. A. Terrero-Escalante and A. A. Garcia, 
Phys. Lett. B517, 243 (2001), astro-ph/0106020. 
I. S. Gradshteyn and I. M. Ryzhik, Tables of Integrals, 
Series and Products, Academic, New York, 1981. 
CETAMA (Commission d'etablissement des methodes 
d'analyses du commissariat a l'energie atomique), Statis- 
tique appliquee a I 'exploitation des mesures" , Masson, 
Paris, 1986. 

S. Leach, A. R. Liddle, J. Martin and D. J. Schwarz, 

Phys. Rev. 66, 023515 (2002), astro-ph/0202094. 

C. P. Burgess, J. M. Cline, F. Lemieux and R. Holman, 

JHEP 0302, 048 (2003), hep-th/0210233. 

N. Kaloper and M. Kaplinghat, Phys. Rev. D68, 123522 

(2003), hep-th/0307016. 

J. Martin and P. Peter, Phys. Rev. D68, 103517 (2003), 
hep-th/0307077. 



